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Abstract 

For uniformly chosen random a G [0, 1], it is known the probability 
the n*'^ digit of the continued- fraction expansion, [a]„ converges to 
the Gauss-Kuzmin distribution P([a]„ = k) ^ log2(l + l/k{k + 2)) 
as n ^ oo. In this paper, we show the continued fraction digits of 
y/d, which are eventually periodic, also converge to the Gauss-Kuzmin 
distribution as d — > oo with bounded class number, h{d). The proof 
uses properties of the geodesic flow in the unit tangent bundle of the 
modular surface, T^{SL21j\M). 

1 Continued Fractions... 

For any a £ [0, 1] we can define the continued fraction expansion in 
Z by repeating a two step algorithm. First ao = a and 6o = [oi\ . Now 
we simply repeat: 

afc+i = {l/fflfe} bk+i = [ak\ (1) 

The end result is that a can be encoded as a sequence of integers: 
[60,^1,^2, •••]• 

If a is rational then we get a finite continued fraction. What if a 
is the square root of a irrational number? Then we get an eventually 
repeating sequences of numbers bk- For example, the sequence for V7 
is [2, 1, 1, 1,4, 1, 1, 1,4, ... ] where the [1, 1, 1,4] motif repeats forever. 
How can we get a purely periodic sequence? A theorem by Galois 
says: 
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Theorem 1.1 (5). A quadratic number a has a purely periodic con- 
tinued fraction expansion if and only if a > 1 and — 1 < a' < where 
a, a' have the same quadratic equation. 



Now let's ask about statistics of these continued fractions. How 
often does the number 5 appear in a generic continued fraction? This 
answer for a random a G [0, 1] chosen uniformly was found by Kuzmin 
in 1928. He showed: 

Theorem 1.2 (5). There exist positive constants A,B such that 
Where An{k) = \{xe [0, 1] : bn{x) = k}\. 



< 



A 



k{k + l) 



In this paper we look at how the statitics of the continued fraction 
digits of ^/d for d > behave as d ^ oo. In fact we have to be more 
specific and restrict ourselves to the case of bounded class number, so 
h{d) is less than some constant. Also, since our sequence 6fc(\/d) is 
deterministic, we need to define the statistics we'll be looking at: 

#{0 < i < T : bija) = k} 
da, k) = lim — — — 

We claim that these statistics approach the limit above, i.e. 
Theorem 1.3. As d — > cxd with h{d) bounded: 

Jimc(v^,fc)^log, + 

To prove this we're going, as d ^ oo and h{d) = 1, the orbits of ^/d 
under the map T : x {1/x} approach the Gauss-Kuzmin on [0, 1]. 
We can rephrase Theorem 1.4 in this new language: 

Theorem 1.4. Let a:;o = {Vd}, h{d) = 1, T : x ^ {V^} be the 
Gauss map and / : [0, 1] — M be continuous: 

fix) dx 



hm hm 1 f(T'{xo)) = [ 

d— *oo AT— »oo iV ^ — ' In 

k=0 

where h{d) = 1 as d goes to infinity. 



- l + x 
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To prove this we need to change settings and examine geodesies in 
the upper half plane. 



Let's switch contexts to H = {x + iy : y > 0} as a differentiable 
manifold with the Poincare metric: 



Much of this exposition follows [2], Chapters 3 and 13. 

The geodesies in this metric arc (Euclidean) semi-circles with diam- 
eters along the real line. Thus for any unit tangent vector {z,e'^^) G 
r^(lHI) there is a unique oriented geodesic which goes through z and 
whose tangent at z points in the direction e*^. 

The group SL(2,M) acts on H by fractional linear transformations: 



Then points in the complex plane are identified as points in the pro- 
jective complex line P-'^(C: 



This is simply the usual matrix action of PSL(2,]R). 

We can also consider the quotient group of H under the action 
of SL2(Z). The quotient under this group action SL2(Z)\]HI has is 
the fundamental domain represented by the intersection of four sets 
{Im(z) > 0}, {\z\ < 1}, {\z + 1| < 1} and {\z - 1| < 1}. 

The tangent space of H is simply H X C. Any element g G SL(2,M) 
can act on the tangent bundle by: 



2 ... and the Geodesic Flow 
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It turns out this action is simply transitive and therefore 

T^{m) ~ PSL(2,R) 
where (H) is the unit tangent bundle. 

Now any two points, zi,Z2 G H determine a unique geodesic - the 
unique Euclidean semi-circle passing through zi and Z2- We can con- 
sider a map, Qt, which flows a tangent vector along its geodesic exactly 
t units of arc length. This is the geodesic flow from zi to Z2- You need 
to specify both a starting point an a direction, an element of S^, to 
get a unique geodesic. Equivalently we can describe the geodesic as 
right-multiplication by the elements: 



This will be the basis for defining the continued fraction map in 
terms of the geodesic flow. Furthermore we can define the geodesic 
flow restricted to T^(SL2(Z)). 

In light of Theorem 1.1, we should consider geodesies whose end- 



— 1 < a' < 0. These curves necessarily cut the set {yi : < y < 1} 
transversely. In fact, we can identify these geodesies either by their 
endpoints or by the tangent vector at which they cut [0, i]. 

The elements oi B = T^[{yi : < y < 1}) parameterize these 
geodesies by the angle at which they cut the line segment [0, i]. If we 
only consider the first case, the set is called S+ and in the second case 
it is called B~. 



These correspond to purely periodic continued fractions and therefore 
to closed geodesies in the Riemann surface SL2(Z)\1HI. 

Definition 2.1. Consider a geodesic, 7 in H with endpoints a, a' for 
which a < — 1 and 1 > a' > 0. Define natural coordinates by {y,z) 
where y = a and z = , . 




e-*/2 
e*/ 





some d > and such that a > 1 and 
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Definition 2.2. Define T : B ^ B in terms of the geodesic flow Qt 
by: 

T[{^, e'^)] = Gtoiz, e'^) where to = M{t > : Qtiz, e'^) e SL2(Z)(S)} 

This to may not be finite but whenever it is finite this map is well- 
defined. This is known as the return time map for the cross section 
B. 

Lemma 2.1. Let x = (16, e*^) G have natural coordinates {y,z). 
Then T{x) G SL2(Z)(5) if T is defined on x. Moreover T{x) G 
SL2(Z)(x'), where x' has natural coordinates 

T{y,z)= 

A similar property holds for 



The first step in the proof of Theorem L3 is to consider closed 

geodesies in SL2(Z)\EI of discrimant d. They correspond to purely 
periodic continued fractions and to elements of B. Their natural co- 
ordinates lie in the ring Q[\/d]. Thus, we can prove these geodesies 
become equidistributed in B as d — oo and h{d) < M, then the first 
natural coordinate follows the Gauss-Kuzmin distribution in [0, 1] . In 
other words, the orbit under T : x ^ {1/x} of any element of Q['\/d] 
becomes Gauss-Kuzmin in [0, 1], asymptotically 



3 Duke's Theorem 

We to define some special sets of geodesies: 

Definition 3.1. Let d < and let (x^, Vd) be the fundamental solution 
to — dy^ = 4. Define Td as the set of geodesies in PSL2(Z)\H of 
length d induced by quadratic forms q{x, y) = ax^ + bxy + cy^ with 
6^ — 4ac = d. 



We can consider closed geodesies on SL2(Z)\H 

Theorem 3.1 (1). Suppose is a fundamental discriminant. Then 
for some 6 > depending only on ft 

1^1 = f^i^) + 0{d-') (2) 
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as d — oo whcc \C\ is the non-Euclidean length of C and the O 
constant depends only on 6 and il. 

The idea that geodesic orbits in homogeneous spaces become equidis- 
tributed can be extended to the tangent bundle. In this case, we con- 
sider the geodesic flow in T^SL2(Z). This is useful because the set we 
wish to consider B,B^ and , who live in the unit tangent bundle 
and not the underlying space. Fortunately for us, Duke's theorem 
extends to this case well. 

Definition 3.2. Define = {j^g^ : [q] e PSL2(Z)\Qd(^)}, the 
geodesies associated with the set of binary quadratic forms modulo 
PSL2 equivalence. If q{x,y) = ax^ + hxy + cy^, the geodesic 7[g] has 
endpoints defined by g(x, 1) = 0. Then project this geodesic onto the 
modular surface PSL2(Z)\H. 

Theorem 3.2 (4). As d 00, d = 0, 1 mod 4, d not a perfect 
square, the set r,^ becomes equidistributed on the unit tangent bundle, 
r^SL2(Z), with respect to the volume measure cZ/xl = 

Now we are ready to prove our equidistribution result. Theorem 1.4 
can be rephrased in dynamical into ergodic theory language. Proving 
it for f{x) = xi{x) for some interval I G [0, 1], we can prove it for any 
continuous function /(x). 

Theorem 3.3. Let T : [0, 1] [0, 1] be the map defined in Theorem 
2.1. Then the orbit {T'"(x) : n G N} is distributed as the Gauss- 
Kuzmin distribution as d goes to infinity. Specifically, for any interval 
/C [0,1]: 

lim - " < • ^ _ 1 f 



d^oo l{d) ln2 1 + X 

Where l{d) is the period of the continued fraction corresponding to 
Vd. 

Proof. Let / C [0, 1] be an interval and (y, z) be the natural coordi- 
nates in B. 

Bi,, = {{y,z) e B : y e 1} X [-e/2,e/2] 
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This set can be embedded in T^{SL2{Z)\M) as: 

{gt{B):te[-e/2,e/2]} 

The Haar measure in the natural coodinates is Lebesgue in all three 
variables. 

By Duke's theorem, there is (5 > such that 

heroic n D,,\ _ 



fi{Bj^,) + 0{d-') (3) 



Note that this equation is true even if / = [0, 1] . It therefore follows 
that: 



V irnB \-u(n ^+^(^ ) (4) 

Because of how we defined -B/^g, the total length is simply the number 
of times the geodesies cut I times e: 

= e • #{0 < n < l{d) : r"(xo) G 1} 

Here l{d) denotes the period of the continued fraction with respect to 
l{d). In the case / = [0, 1] the last line is just e • l{d) so the left hand 
side of (4) is really just counting measure: 

#{0 < n < l{d) : T"(a;o) G 1} 

m 

For the right hand side of (4) let's find the measure of Bi^^: 

Therefore equation (4) should read: 

#{0 < n < l{d) : r"(aro) e i"} _ 1 f dy 



lim 

(i-»oo l{d) ln2 J J 1 + y 

□ 



ln2y, 1 



That concludes the proof that these continued fractions of quadratic 
irrationalities follow the Gauss-Kuzmin distribution as the discrimi- 
nant tends to infinity. Next we show this results describes truly generic 
behavior. 
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4 Bounded Class Number 



We can relax the condition h{d) = 1 in Theorem 1.4: 

Theorem 4.1. Let xq = {Vd}, h{d) = I, T : x ^ {l/x] be the 
Gauss map and / : [0, 1] ^ M be continuous: 

^ f{x) dx 



lim lim y^f{T\x^))= / 



d-^oo n— »oo In. In 2 1 + X 

fc=0 

where there exists M > such that h{d) < M. 

Proof. We simply examine Theorem 3.3 more closely. In this case, T^ 
has more than one element, i.e. there are several geodesies with the 
same discrimant d. Duke's theorem says the sum 

yZcpA \cnn\ yZr^A \cnn\ , 

is the Lebesgue measure. Here is the length of a geodesic of dis- 
criminant d. Following the principle outlined in [6] (Section 1.3.5 (1)), 
since the Lebesgue measure is an extreme point in the convex space of 
geodesic-flow invariant measures on T^(SL2(Z)), since the Gauss map 
T is ergodic and since their sum is a Lebesgue measure, each term in 
the sum must also approach Lebesgue measure. □ 



Then continued fractions in Q[vd] tend towards the Gauss Kuzmin 
distribution even in the case bounded class number. This behavior is 
truly generic. 
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